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ABSTRACT 

We consider the production of 3 He and 2 H by 4 He photodisintegration initiated by non- 
thermal energy releases during early cosmic epochs. We find that this process cannot be 
the predominant source of primordial 2 H since it would result in anomalously high 3 He/D 
ratios in conflict with standard chemical evolution assumptions. We apply this fact to 
constrain topological defect models of highest energy cosmic ray (HECR) production. Such 
models have been proposed as possible sources of ultrahigh energy particles and 7-rays with 
energies above 10 20 eV. The constraints on these models derived from 4 He-photodisintegra- 
tion are compared to corresponding limits from spectral distortions of the cosmic microwave 
background radiation (CMBR) and from the observed diffuse 7-ray background. It is shown 
that for reasonable primary particle injection spectra superconducting cosmic strings, unlike 
ordinary strings or annihilating monopoles, cannot produce the HECR flux at the present 
epoch without violating at least the 4 He-photodisintegration bound. The constraint from the 
diffuse 7-ray background rules out the dominant production of HECR by the decay of Grand 
Unification particles in models with cosmological evolution assuming standard fragmentation 
functions. Constraints on massive black hole induced photodisintegration are also discussed. 



1 Introduction 



In this paper we consider various constraints inferred from the possible photodisintegra- 
tionof 4 He in the early universe. Following Protheroe, Stanev, and Berezinsky M we note 
that the photodisintegration of this isotope can be employed to place stringent limits on early 
cosmic energy injections associated with, for example, decaying particles p], |§, evaporating 
black holes [[|, or annihilating topological defects || |6], |7|, |8|, |9], [L0|. Our focus here will 
be particularly on constraining the latter scenario. It has also been suggested that 4 He- 
photodisintegration in the early universe could be a production mechanism for the observed 
light-element abundances of deuterium and 3 He ||11|| . In this work we will study the feasibility 
of such a scenario and show that the ( 3 He/ 2 H) ratio poses a problem to it. We will show 
that photodisintegration yields ( 3 He/ 2 H) >> 1 and since 2 H is destroyed and 3 He increases 
with evolution, measures of ( 3 He/ 2 H) place severe constraints on photodisintegration. 

Nonthermal energy releases at high redshifts may leave various observable signatures. 
The cosmic microwave background radiation (hereafter, CMBR) has been measured to have 
a blackbody spectrum to very high accuracy [|12[]. Any injection of energy between redshifts 
of z ~ 10 3 and z — 3 x 10 6 may produce observable spectral distortions of the blackbody spec- 
trum Here the lower redshift represents the approximate epoch of decoupling (assuming 
no re-ionization), whereas the higher redshift represents the epoch at which double-Compton 
scattering is still efficient enough to completely thermalize significant energy releases |TJ] . 



The diffuse 7-ray background observed at the present epoch can also be used to constrain 
early cosmic energy injections [PJ. For redshifts z < 300 — 1000 pair production by 7-rays 



on protons and 4 He is rare so that the universe becomes transparent to 7-rays with energies 
below -Emax- Here the energy E max is 

Bn ^_^ 17GeV (_) , (!) 

where T is the CMBR temperature and m e is the electron mass. E max is related to the 
threshold energy for e + e~-pair creation by high-energy 7-rays scattering off CMBR-photons. 
Any radiation with energies above this threshold is effectively instantaneously "recycled" by 
pair production (77cmbr — ► e + e~) and inverse Compton scattering of the created electrons 
and positrons (c7cmbr — ► 67). These processes yield a degraded 7-ray spectrum with generic 
energy dependence cx E~ 15 considerably below E mSuX before steepening and finally cutting off 
at E max 0. Significant energy releases in form of high-energy 7-rays and charged particles 
at epochs with redshifts below z ~ 300 — 1000 may therefore produce a present day 7-ray 
background and are subject to constraint. 

For redshifts smaller than z — 10 6 stringent constraints on various forms of injected 
energy can also be derived from the possible photodisintegration of 4 He and the concomitant 
production of deuterium and 3 He. The injection of high-energy particles and 7-rays above 
the energy threshold E max will initiate an epoch of cascade nucleosynthesis subsequent to 
the epoch of standard primordial nucleosynthesis at T ~ 100 keV. The abundance yields of 
2 H and 3 He produced by 4 He-photodisintegration during cascade nucleosynthesis are quite 
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independent from the primary 7-ray and charged particle energy spectra. Deuterium and 
3 He abundance yields depend only on the amount of injected energy and the injection epoch. 
For the detailed calculations leading to these conclusions the reader is referred to the work 
by Protheroe, Stanev, and Berezinsky The nucleosynthesis limits on the release of energy 
into the primordial gas can be up to a factor of ~ 100 more stringent than equivalent limits 
on energy releases derived from distortions of the CMBR-blackbody spectrum. 

For redshifts z > 10 6 , corresponding to CMBR-temperatures of T > 200 eV, the photo- 
disintegration of 4 He is inefficient. This is because the energy threshold for pair production 
falls below the energy threshold for 4 He-photodisintegration , -E max < E t ^ e . The best nu- 
cleosynthesis limits on decaying particles and annihilating topological defects in the cosmic 
temperature range IkeV < T < lOkeV come from the possible photodisintegration of deu- 
terium |3], [L6| . These limits are stronger than analogous limits from distortions of the CMBR 
blackbody spectrum. 

In this narrow temperature range limits on decaying particles and topological defects 
may, in fact, be more stringent due to effects of injecting antinucleons. Antinucleons may 
be produced during 77cmbr pair production for 7-energies > 10 5 GeV or when there is a 
significant hadronic decay channel for a massive decaying particle or topological defect. These 
antinucleons can then annihilate on 4 He thereby producing approximately equal amounts of 
2 H and 3 He ||17||. We will, however, not further pursue this idea here. 



For temperatures above T ~ 1 keV there are virtually no constraints on decaying parti- 
cles and topological defects from distortions of the CMBR blackbody spectrum. However, 
stringent limits on decaying particles and topological defects may obtain from the injection 
of hadrons (for a review see |9). An injection of mesons and baryons generally increases the 
neutron-to-proton ratio and results in increased 4 He-mass fractions (1 MeV > T > 100 keV) 
and/or increased 2 H and 3 He- abundances (100 keV > T > lOkeV; [0). It has been sug- 
gested that a combination of 4 He-hadrodestruction and 2 H, 3 He-photodestruction induced 
by a late-decaying particle (T ~ 3 keV) may bring big-bang-produced light-element abun- 
dances close to observationally inferred abundance constraints for a wide range of fractional 
contributions of baryons to the closure density, Qt, |1| . 

The observational signatures of such scenarios are primordial isotope ratios of ( 3 He/ 2 H)~ 
2 — 3 and 6 Li/ 7 Li ~ 1, contrasting the predictions of a standard, or inhomogeneous, big-bang 
freeze-out from nuclear statistical equilibrium. For a wide range of parameters, such as de- 
caying particle life times and hadronic branching ratios, these models would overproduce 2 H 
and 3 He and therefore the calculations by Dimopoulos et al. [19] do also serve as constraints 
on particle parameters and abundances. We note here that the high ( 3 He/ 2 H) ratio may in 
fact be a severe problem for such scenarios. 

In this paper we restrict ourselves to constraints derived from the effects of nonthermal 
energy injections at epochs with redshifts z < 10 6 . The outline of the paper is as follows. 
In Section 2 we briefly review the observationally inferred light-element abundances of 2 H 
and 3 He. We then consider 4 He-photodisintegration scenarios and their compatibility with 
the observations. In Section 3 we study the effects of possible energy injection by super- 
conducting strings, ordinary strings, and magnetic monopoles on the primordial 2 H and 3 He 
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abundances, the distortions of the CMBR-blackbody, and the diffuse 7-ray background. In 
these scenarios we assume that such topological defects would radiate on a level such that 
they could produce the observed highest energy cosmic rays at the present epoch. Conclu- 
sions are drawn in Section 4. Throughout this paper we will mostly use c = h = 1. 

2 Constraints on 4 He-Photodisintegration as the pre- 
dominant Source of Primordial Deuterium 

In this section we investigate scenarios which have 4 He-photodisintegration as an efficient 
production mechanism of the light-element abundances of deuterium and 3 He. In this study 
we are naturally led to consider the primordial ratio of ( 3 He/ 2 H) p . This is because the ratio 
of these light isotopes emerging from the big bang nucleosynthesis process, ( 3 He/ 2 H)BBN , 
is quite different from that emerging from the 4 He-photodisintegration , ( 3 He/ 2 H) p hoto- In 
particular, we expect generic isotope ratios of ( 3 He/ 2 H)BBN ^ 1, and ( 3 He/ 2 H) p hoto 3> 1- 
We will show that this fact can be used to severely constrain the photodisintegration of 4 He 
as the principal source of primordial deuterium. We will also show that the observationally 
inferred abundances of 2 H and 3 He may imply a factor 2-3 more stringent constraints on the 
primordial number densities of decaying particles and on the energy injected by topological 
defects than previous work has assumed. 

The most accurate determination of a ( 3 He/ 2 H)-ratio is thought to come from solar 
system observations of 3 He abundances. Geiss |OT reanalyzed the existing data and inferred 
for the abundances of deuterium and 3 He at the time of solar system formation 

1.2xl0- 5 < C-^-) < 1.8 xlO- 5 , 



H Jr. 







1.6xl0- 5 < (y) £ 3.3xl0~ 5 , (2) 
' 3 He\ 

0.34 < [— < 1.13 . 







A determination of the interstellar medium abundances of 2 H and 3 He is less precise due 
to observational difficulties [2TJ. The observed ( 2 H/H)-ratios ranges between 5 x 10~ 6 < 



( 2 H/H)ism ^ 2 x 10 5 [p2| . Interstellar ( 3 He/H)-ratios are observed in the range 1.1 x 10 5 < 
( 3 He/n)isM < 4.5 x 10~ 5 [p3"| . These abundances imply a present ( 3 He/ 2 H)-isotope ratio of 
0.55 < ( 3 He/ 2 H) ISM < 9. 

Deuterium is the most fragile of the light isotopes. It is easily destroyed during the 
pre-main sequence evolutionary stage of stars via 2 H(p,7) 3 He. Furthermore, there are no 



plausible galactic production sites for deuterium. Epstein, Lattimer, and Schramm p4 
summarize the arguments against a galactic origin of deuterium. The chemical evolution 
of 3 He is less clear. It is known that 3 He is destroyed to some extent in massive stars 
(M > 5 — 8M ), whereas low-mass stars (M < 1 — 2M Q ) may be net producers of 3 He. 
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This theory is supported by the observations of 3 He abundances in planetary nebulae. It 
is certainly very reasonable to assume that standard chemical evolution models can only 
increase the primordial ( 3 He/ 2 H) p -ratio, 

~He\ / 3 He\ 



> 

2 H I ■ ~ ^ 2 



In this expression ( 3 He/ 2 H) t denotes the isotope ratio at some cosmic time t and the pri- 
mordial isotope ratio ( 3 He/ 2 H) p includes any pre-galactic production mechanism, such as 
big bang nucleosynthesis and 4 He-photodisintegration in the early universe. Note that the 
inferred ( 3 He/ 2 H) ratios at the time of solar system formation and the present epoch are 
consistent with the assumption of monotonically increasing ( 3 He/ 2 H) ratios with time. 

The ( 3 He/ 2 H)-ratio in a standard homogeneous big bang nucleosynthesis (hereafter, 
SBBN) scenario at baryon-to-photon ratio 77 = 3 x 1CT 10 is ( 3 He/ 2 H)sBBN — 0.2. An upper 
limit on the ( 3 He/ 2 H)-ratio in SBBN can be obtained by requiring the 4 He-mass fraction 
to satisfy Y p < 0.25, whereas a lower limit on this isotope ratio can be estimated from the 
conservative bound ( 2 H/H) < 3 x 10 -4 . This yields the SBBN range 

0-09 < (5^) < 0.55. (4) 

V Z M / SBBN 

Typical ( 3 He/ 2 H)-isotope ratios resulting in inhomogeneous big bang scenarios are not very 
different from those in Eq. 

The detailed calculations by Protheroe, Stanev, and Berezinsky show that the abun- 
dance ratios of ( 3 He/ 2 H) produced during cascade nucleosynthesis in the early universe 
exceed 

W) ^ 8 > (5) 

tl / photo 

for a wide range of fractional contributions of baryons to the closure density, Hubble 
parameters Hq in units of 100 km sec -1 Mpc -1 , h, and epochs of energy injection. This is 
because in 4 He-photodisintegration the effective cross sections for the two-nucleon photoab- 
sorption processes [ 4 He(7,pn) 2 H and 4 He(7, 2 H) 2 H] are roughly ten times smaller than the 
effective cross sections for the single-nucleon photoabsorption processes [ 4 He(7,p) 3 H and 
4 He(7,n) 3 He] ||. 

Note that Eq. (|5|) applies strictly only under the following assumption. In cascade nu- 
cleosynthesis it is assumed that the main fraction of radiation is injected above the energy 
threshold Eq. ([!]) for 77cmbr — > e~e + pair creation. Pair creation and inverse Compton 
scattering will then yield a generic 7-ray spectrum with energy dependence oc _E~ L5 be- 
low E max /2 and oc E~ b above before cutting off at E ma _ x . These 7-rays can be effective 
in photodisintegrating 4 He where the competing process is the consumption of 7-rays by 
Bethe-Heitler pair production on hydrogen and helium. 

When radiation is injected below -E max the 7-rays may have a spectrum quite different 
from the behavior oc E~ l,b depending on the actual 7-ray source. In principle, it is then 
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conceivable to photodisintegrate 4 He in such a way that isotope ratios of ( 3 He/ 2 H)~ 1 result. 
This could be accomplished by a 7-ray source which preferentially radiates above energies of 
E ~ 100 MeV but below E max . This is because only in the energy range between the 4 He- 
photodisintegration threshold E t ^ e = 19.8 MeV and E ~ 100 MeV the effective cross section 
for 3 He production in 4 He-photodisintegration is roughly ten times larger than the effective 
cross section for 2 H production in this process. For 7-ray energies E > lOOMeV these cross 
sections are roughly equal. In practice, any such scenario has to occur at relatively low 
redshifts z < 10 3 so that there will not develop a "softer" second generation 7-ray spectrum 
produced by Bethe-Heitler pair production and inverse Compton scattering. In this case, 
however, significant deuterium production would require 7-ray fluxes which would exceed 
the present day diffuse 7-ray background. 

It should be noted that 7-rays could also be effective in photodisintegrating 3 He and 
2 H and thereby in resetting any initial ( 3 He/ 2 H) photo -isotope ratio produced during cascade 
nucleosynthesis to smaller values. However, the relative abundances of 4 He-targets to Re- 
targets is approximately 10 3 — 10 4 to 1, so that for roughly equal photodisintegration cross 
sections the number densities of 7-rays in the energy range between the 3 He-photodisintegra- 
tion threshold E^ e = 5.4 MeV and E^ e = 19.8 MeV should be 10 3 - 10 4 times larger than 
the number densities of 7-rays with energies above E t ^ e . Such a scenario would require an 
extremely "soft" 7-ray spectrum. 

We can derive limits on the allowed contributions of 4 He-photodisintegration to the pri- 
mordial 2 H and 3 He abundances. This can be done by employing the solar system ( 3 He/ 2 H)- 
isotope ratio from Eq. (§) and assuming that this ratio represents a conservative upper limit 
on the primordial ( 3 He/ 2 H)-isotope ratio [refer to Eq. (|^)]. Note that when either one of 
Eqs. (0) or @ does not apply one of the widely used standard assumptions of galactic chem- 
ical evolution has to break down. We can derive an upper limit on the fraction of deuterium 
/ 2 P H° to contributed to the primordial deuterium abundance by 4 He-photodisintegration. A 
simple calculation of the abundance average then yields 

/ 3 He \ _ / 3 He \ 
/•photo < - I 2 H )q \ 2 H JBBN 

J2 H ~ /£Hex _ /3 He \ • W 

I 2 H Jphoto I 2 H JbBN 

By using the upper limit for ( 3 He/ 2 H) Q from Eq. the lower limit in Eq. (|J) for the 
( 3 He/ 2 H) SB BN- r atio, and Eq. (|5|) for the ( 3 He/ 2 H) photo -ratio we derive 

^photo < 13% _ (7) 

It is evident that the contribution of deuterium produced during cascade nucleosynthesis to 
the total primordial deuterium abundance has to be small in order to not overproduce 3 He. 
This also implies that generic 4 He-photo disintegration scenarios can not be the predominant 
production mechanism of the primordial 2 H and 3 He light-element abundances. The stringent 
limit of Eq. (|7|) can only be evaded when either there existed an extremely "soft" 7-ray source 
in the early universe or when generic features of the galactic destruction/production of 3 He 
and 2 H are for some yet unknown reason not understood. 



5 



Gnedin and Ostriker JTT[] have proposed the interesting scenario of a very early formation 
(z ~ 800) of massive black holes. If these black holes do accret material which emits a 
quasar-like X-ray and 7-ray spectrum they may induce the photodisintegration of 4 He and 
the reionization of the universe. The reionization of the universe would cause primordial 
CMBR fluctuations to be erased, whereas the processed 7-ray spectrum could constitute 
the diffuse 7-ray background at the present epoch. They concluded that this selfconsistent 
model could evade the upper limit on Q b given by the observed deuterium abundance and a 
SBBN scenario since deuterium and 3 He would have been produced, at least in part, in the 
4 He-photodisintegration process. For typical models they produce a fraction i?^ * — 50% of 
the total primordial deuterium abundance by 4 He-photodisintegration. Clearly, this fraction 
is in conflict with the limit of Eq. ([?]) and would result in too high ( 3 He/ 2 H) p ratios [ 26| . 



fft%^ < 35% - 55% . 



We can also constrain the fraction f?2^3 Be \ which can be contributed to the total sum of 
the primordial deuterium- and 3 He-abundances by 4 He-photodisintegration. This parameter 
is limited by 

'( 2 H+3He) 

Any annihilating topological defects or decaying particles abundant enough to initiate an 
epoch of cascade nucleosynthesis such that more than 35% of the presently observed abun- 
dance sum of ( 2 H+ 3 He) is contributed by this cascade nucleosynthesis are subject to con- 
straint. The limits given in Eq. (|8|) are a factor 2-3 better than equivalent limits assumed in 
previous work. 

These limits can be put into context by the upper limit on the sum of 2 H and 3 He inferred 



from the solar system data and chemical evolution models by Geiss |20 

- 2 H+ 3 He 



H 



<l.lxl0- 4 . (9) 



In Figure 1 we show constraints from 4 He-photodisintegration on the maximum allowed en- 
ergy release as a function of redshift. To produce this figure we have used Eq. (|9|) and 
the upper range given in Eq. (||). For comparison we show analogous limits from possible 
distortions of the CMBR background. These are taken from reference El. It is seen that 



over a wide range of redshifts the limits from 4 He photodisintegration are more stringent 
than the limits from CMBR distortions. Also shown are constraints from the diffuse 7-ray 
background which result from the generic cascade spectrum (see section 3.4). 

3 Energy Injection from Topological Defects 
and Highest Energy Cosmic Rays 

3.1 History of Energy Injection in Defect Models 

It is commonly believed that cosmic rays are produced mostly by first order Fermi accel- 
eration (see e.g. |27|, |28|) at astrophysical shocks in the presence of magnetic fields. The 
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highest energies seem to be reached in relativistic shocks contained in radiogalaxies and ac- 



tive galactic nuclei (see e.g. |29|, |3T], |32|j). The recent observation of cosmic rays above 
10 20 eV by the Fly's Eye |33], EM and AGASA GO, ETH] experiments, and the experiment 



at Yakutsk |38|, |39| may, however, not be easily explained by this mechanism [^, [41, [42 



Therefore, it has been suggested that such superhigh energetic cosmic rays could have a 
non-acceleration origin |5|, [|, |41], [E| [44], ^5], |46| as, for example, the decay of supermassive 



elementary "X" particles associated with Gand Unified Theories (GUTs). These particles 
could be radiated from topological defects (TDs) formed in the early universe during phase 
transitions caused by spontaneous breaking of symmetries implemented in these GUTs. This 
is because TDs, like ordinary or superconducting cosmic strings and magnetic monopoles, 
on which we will focus in this paper, are topologically stable but nevertheless can release 
part of their energy in form of these X-particles due to physical processes like string collapse 
or monopole annihilation. The X-particles with typical GUT scale masses (~ 10 15 GeV) 
decay subsequently into leptons and quarks. The strongly interacting quarks fragment into 
a jet of hadrons which results in typically of the order of 10 4 — 10 5 mesons and baryons. 
It is assumed that these hadrons then give rise to a substantial fraction of the HECR flux, 
whereas the contribution from the lepton primary is often approximated to be negligible. 
It also causes a more or less uniform global energy injection whose spectrum is determined 
by the cascades produced by the interactions of the primary decay products with various 
background radiation fields. This energy injection is subject to the constraints from 4 He- 
photodisintegration discussed in the previous section as well as to constraints from spectral 
CMBR distortions and the observed 7-ray background. 

The X-particle injection rate dnx/dt as a function of time t or redshift z usually is 
parametrized as |44| 

^oct-^. (10) 
at 

It is important to note that the effective value of p may depend on the epoch. Given that and 
using standard cosmological relations for t(z) [47] one can describe the X-particle injection 



history by introducing the dimensionless function 

^ (dn x /dz)(z) 
t (dn x /dt)(t ) ' 

where to = 2Hq X /3 is the age of the universe (we assume a flat universe, f2 = 1, throughout 
this paper). 

For example, for annihilating magnetic monopoles it can be shown that p = 1 for 
t > t eq and p = 1.5 for t < t eq , where t eq is the time of matter-radiation equality. 

As a second example, let us look at collapsing cosmic string loops. These may include 
ordinary as well as superconducting strings. Let us assume that the history of loops consists 
of two distinct evolutionary stages. We will see below that such a schematic representation 
can be used for both superconducting strings and ordinary strings. In the first stage the 
loop slowly radiates gravitational radiation with a power ~ lOOG/i 2 . Here, G is Newton's 
constant and fi ~ v 2 is the energy per unit length of the string in terms of the GUT symmetry 
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breaking scale v. This will decrease the loop length, Lit), at an effective rate v g ~ lOOG/z, 



L(t) = L b -v g (t- t b 



(12) 



In this expression tb and L& denote the birth time and the loop length at birth, respectively. 



Numerical string simulations |I8, 49], [5(J suggest that loops are born with a typical length 



Lb = atb with a being a dimensionless constant which can be as small as a few times v g . To 
simplify the calculation we will assume that all loops are born with the same length Lb. 

Note that the gravitational radiation associated with this first stage of string loop evo- 
lution should not have any effects on CMBR distortions, the diffuse 7-ray background, or 
result in photodissociation of 4 He. The existence of gravitational radiation during the epoch 
of primordial nucleosynthesis, however, can effect abundance yields by changing the cosmic 
expansion rate flSlf . For symmetry breaking scales v < 10 16 GeV as discussed in this paper 
this effect is negligible. 

Once the loop enters the second evolutionary stage gravitational radiation becomes a 
subdominant energy loss mechanism. The loop starts to collapse at a rate which grows 
considerably beyond the gravitational rate v g by radiating other forms of energy, one of 
them being X-particles. The decay products of these X-particles may then contribute to the 
HECR flux observed at the present epoch. We schematically assume here that during this 
second evolutionary phase a fraction / of the total energy in loops smaller than a certain 
critical length scale, L c (t), is instantaneously released in form of X-particles. This is a good 
approximation as long as the time which loops spend in their second evolutionary phase is 
short compared to the cosmic time t. Denoting the birth rate of closed string loops per unit 
volume being chopped off of the string network at birth time tb by (drib/dt) tb we can then 
write the rate of X-particle production per unit volume as 



dn x , , _ dn b 
~df [) ~ f ~oT 



dtb 
~dt 



Rjtb) 
R(t) 



yL c (t) 
m x 



(13) 



Here R(t) is the cosmic scale factor and mx = gv is the X-particle mass in terms of the 
symmetry breaking scale v and the Yukawa-coupling g (g < 1). Furthermore, (dtb/dt) takes 
account of the time delay between the birth of a string loop at time tb and the final phase of 
X-particle evaporation at later time t. Finally, the factor [R(tb) / R(t)} 3 accounts for dilution 
due to the cosmic expansion between tb and t. If the string network exhibits scaling behavior 
the birth rate of closed string loops can be written as 



dnb 
~~dt 



1 

tV 



(14) 



where (3 is a dimensionless constant which is approximately related to a by the relation 
a(3 ~ 0.1 P. 

The possible existence of superconducting cosmic strings within certain GUTs was first 
proposed by Witten ||. Ostriker, Thomson and Witten [0 (hereafter, OTW) discussed quite 
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severe potential cosmological consequences and also suggested that these objects might con- 
tribute to the ultrahigh energy cosmic ray flux. This was further pursued by Hill, Schramm, 
and Walker || who mainly investigated fermionic superconducting string loops which could 
produce HECR by ejecting superheavy fermion pairs towards the end of their evolution. 
With respect to the schematic scenario described above two cosmic epochs have then to 
be considered for superconducting cosmic strings of this type. For cosmic time t < t tr , all 
existing loops are radiating dominantly in electromagnetic and/or X-particle radiation. In 
terms of the (in general time dependent) saturation length L s (t) the relevant condition which 
compares electromagnetic and gravitational energy loss rates reads f| 



gL s {ty 2 

at 



>v g . (15) 



In this case the loops have not experienced a gravitational radiation dominated energy loss 
phase, but rather have directly entered the phase of comparatively fast collapse at birth. 
We can therefore approximate ~ t and the critical length L c (t) is the minimum of the 
birth length, Lb(t), and the saturation length, L s (t). In contrast, for cosmic times t > t tr the 
epoch at which a string loop reaches its second evolutionary stage is primarily determined 
by gravitational energy loss, U — (v g /a)t, and L c (t) = L s (t). Using Eqs. ( P~3|) and (|nj) this 
leads to the following time dependence of the X-particle injection rate: 



dnx^j r 4 ™p L s {t) if t>t 



R(t) 



dt 1 t- 4 Mm(L s (t),at) if t < t 



(16) 



The saturation length for superconducting strings depends on the intergalactic magnetic 
field history || and is therefore strongly model dependent. OTW originally considered an 
intergalactic field whose energy density scales like the CMBR energy density. In this scenario 
the saturation length is roughly constant in time and can be written as 

L s (t) ~ const. ~ 10 ( -^-] ( -A-V ( " V V Vv/ 3 pc, (17) 
sW V10- 9 G/ V 1M P C / V10 15 GeV/ y P ' v ' 

where B and Ao are strength and coherence length of the intergalactic field today. Using 
Eq (|T5|), the transition time t tT which separate the two string evolution epochs is in terms of 
redshift z tr given by 

Za = 4. 78 x 10 3 ( ^y 1/2 (J^Y f_^Y 2/ V/ 6 



10~ 9 G/ V 1M P C / V10 15 GeV 

For the calculations performed in the following we will use z tr = 2 x 10 3 . Since we will match 
the two functional time dependences in Eq. (|16D at t — t tr and since Min(L s (t), at) = at 
for t ttr we w ih use dn x /dt oc t~ 3 for all t < t tr for a lower bound on energy injection. 
Furthermore, we will neglect the time dependence coming from the factor [R(v g t/a)/R(t)] 3 
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in Eq. flT6|) in case U = v g t/a < t cq and t > t cq . A more detailed treatment would have to 
take into account the chronological order of tb, t and t eq as well as the finite collapse time 
of a string loop in its second evolutionary stage. This would be model dependent via the 
parameters from Eqs. (|T7|) and ([18]). It is, however, easily seen that such effects lead to 
X-particle injection rates which can only be larger at early times than the injection rates of 
our simplified treatment Eq. ([161) . Our calculations will therefore give us conservatively low 
estimates for the total energy release in X-particles. Within these approximations Eq. ( |16|) 
is of the form of Eq. fliCf ) with p = for t > t tT and p = 1 for t < t tT . 

In principle, for superconducting strings the energy radiated in form of X-particles is 
determined by the model. In Ref. |53[ it was shown that the ultrahigh energy particles are 
absorbed in the strong magnetic field produced by the electric current in the string loops. 
Instead, it was suggested that most of the string energy would be liberated in the form 



of neutrinos 54]. We shall demonstrate here that even without these effects in the OTW 



scenario, where L s (t) is approximately constant in time, it is barely possible to produce 
the observed HECR flux for reasonable model parameters. It has been shown || that in 
scenarios where L s (t) grows with time the saturation length at the present epoch, L s (i ), 
has necessarily to be smaller than the L s (to) in the OTW scenario. Such scenarios would, for 
example, be given when intergalactic magnetic fields are increased by dynamo effects. In this 
case it follows from Eq.flTBD that the HECR flux at the present epoch can not be produced 
by superconducting cosmic strings even when / = 1. In the opposite case (L s (t) growing 
in time) scenarios are conceivable where an / < 1 can reproduce the observed HECR flux. 
However, for a given universal HECR flux more energy would have been injected outside of 
the strong magnetic field region in the past compared to the OTW scenario. Therefore, if 
too much energy tends to be injected within the OTW scenario, as will be shown to be the 
case below, the other scenarios are also unlikely to be able to explain the observed HECR 
flux. 

In the case of ordinary strings it has been shown that well known physical processes like 



cusp evaporation are not capable of producing detectable cosmic ray fluxes [II], [55|] . It has, 
however, been suggested that a small fraction / of all loops could be formed in states which 
would lead to their total collapse within one oscillation period after formation The 
total energy in these kinds of loops would be released in form of X-particles. Then, ~ t 
and L c (t) = ~ at so that 



dt 



faP/imxH- 3 , (19) 



which is of the form of Eq. fllUD with p — 1. Recently, there has been a claim |jB that loops 
in high-harmonic states are likely to self-intersect and decay into smaller and smaller loops, 
finally releasing their energy in relativistic particles. Eq. fllpp would be a reasonable good 
approximation also in this case. 

Up to now we have only considered the functional form of the X-particle injection rate 
dnx/dt up to an absolute normalization. If we assume that HECR are produced by decaying 
X-particles radiated from topological defects we can normalize to the differential HECR 
flux j'hecr(-E') observed today (t = t ) at a fixed energy E = E b B . In these models one 
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expects to observe mainly 7-rays at energies E > 10 20 eV [p7|. We define the effective X- 



particle fragmentation function into 7-rays, (dN 1 / dx)(x) where x = 2E/m x , as the effective 
differential primary 7-ray multiplicity per injected X-particle multiplied by 2/m x ||10|| - Then 
the normalization depends on the 7-ray attenuation length X-y(E) and on (dN 1 /dx)(x) at 
x = 2E ohs /m x : 

-1 

jHECR^obs) (20) 



drix 



dt 



(to) 



2iim x 



\ 1 {E t 



obs } 



dN~ (2E ohs 



8.16 x 10 



dx 



-10 



~dN y 


f"2E ohs \ 


dx 


V m x J 



V m x 

10 16 GeV 
1 

cm -3 



A 7 (E obs ) \ 1 / jHECR(-^obs) • GeV cm 2 sec sr 



lOMpc J \ 



4 x 10" 31 



sec 



>obs 



In the last expression of Eq. (|20|) and in the following we have used the numbers for E 
2 x 10 20 eV. 

Using the parametrization of X-particle injection history, Eq. (|ll|), and the normalization 
Eq. (|20|) we are now in a position to derive various constraints on TD models for HECR from 
limits on energy injection into the universe. 



3.2 Limits from Cascade Nucleosynthesis 

In Ref. [|I|] the number iV( 3 He, D, z) of 3 He and D nuclei produced via 4 He-photodisintegra- 
tionper GeV electromagnetic cascade energy injected into the universe was calculated as a 
function of redshift z. These functions depend only weakly on h and fib- Therefore, using 
Eqs. (|TT|) and ( p0|) and assuming that a fraction f c of the total energy release in high energy 
particles goes into the cascade one gets 



Tle\ 
H J 



9.7 f c 



'Q h h 2 



-1 



h 



-1 



photo 



x 0.02 J 

^JHECR(-Eobs) 



\0.75 J 

GeV cm 2 sec sr 



m x 
10 16 GeV ) 



\ 2 ( \\E h$) 



4 x 10" 31 



dN„ (IE, 



dx 



lOMpc 
-1 



v obs 



m x 



AT( 3 He, z) 



l + z) 



(21) 



rdz . 



where the integral is performed over the range in Fig. 4 of Ref. |T|. An analogous formula 
applies for the produced deuterium fraction ( 2 H/H) photo . Using Eq. (§) and the bound 
( 3 He + 2 H)/H < 1.1 x 10~ 4 we can impose the constraint 



' 3 He + D ^ 



< 5 x 10" 



photo 

This leads to lower limits on the fragmentation function taken at x 
the three cases discussed in the previous section read 

1.4 x 10 5 for monopole annihilation 
2.0 x 10 6 for ordinary strings \ x f c 

1.8 x 10 11 for the OTW scenario 



(22) 

2E b s /m x which in 



dNj_ ( 2E ohs 
dx \ m x 



> 



'Q b h 2 
0.02 



h 



-1 



0.75 



11 



m x 



'A 7 (£ D bs)\ 1 /jHECR(^obs) • GeVcm 2 secsr' 



10 16 GeV/ VlOMpc/ \ 4xl0- 31 



(23) 



This has to be compared with expected fragmentation functions in the different defect 
scenarios. In case of monopoles and ordinary strings this function is mainly determined 
by the hadronization of the fundamental quarks created in X-particle decays. At HECR 



energies it is reasonable to assume a power law behavior jig, |46[ . In superconducting string 
scenarios the effective spectrum of HECR, which if at all able to leave the high magnetic field 
region around these strings, could well be altered by interactions with these strong fields. 
Nevertheless it is still reasonable to assume that at least at HECR energies this spectrum 
has a power law form. 

It can easily be shown that a properly normalized power law fragmentation function 
(dN 1 /dx){x) oc x~ q (q > 0) obeys (dN^/dx)(x) < 2x~ 2 for all q > 0. Thus, because of 
Eq. (^3|) the OTW scenario is inconsistent with these power law fragmentation functions 
independent of mx as long as f c > 6.9 x 10~ 3 . In contrast, the monopole annihilation and 
ordinary cosmic string scenarios are compatible with reasonable fragmentation functions. 



3.3 Limits from Cosmic Microwave Background Distortions 

Early non-thermal electromagnetic energy injection can also lead to a distortion of the cos- 
mic microwave background. We focus here on energy injection during the epoch prior to 
recombination. A comprehensive discussion of this subject was recently given in Ref. [|58 |. 
Regarding the character of the resulting spectral CMBR distortions there are basically two 
periods to distinguish: First, in the range 3 x 10 6 ~ z th > z > z y ~ 10 5 between the 
thermalization redshift z t h an d the Comptonization redshift z y , a fractional energy release 
Au/u leads to a pseudo-equilibrium Bose-Einstein spectrum with a chemical potential given 
by yu ~ 0.71 Au/u. This relation is valid for negligible changes in photon number which is a 
good approximation for the Klein-Nishina cascades produced by the GUT particle decays we 



are interested in [58]. Second, in the range z y > z > z rec ~ 10 between z y and the recom- 



bination redshift z rec the resulting spectral distortion is of the Sunyaev-Zel'dovich type EP| 



with a Compton y parameter given by Ay = Au/u. The most recent limits on both /i and y 
were given in Ref. [I2|. The resulting bounds on Au/u for instantaneous energy release as a 



function of injection redshift ]13| are shown as the dashed curve in Fig. 1. 

Since energy injection by topological defects would be a continuous process it is convenient 
to define an effective fractional energy release into the CMBR in the following way: 



Au 

u 



hmx r*dn x t(z) (24) 



off M Jz ICC dz (1 + 2;) 



Here fb is the fraction of the total energy release in high energy particles which contributes 
to the CMBR distortion, u is the CMBR energy density today, and £(z) is given by 10~ 4 
divided by the function shown as the dashed curve in Fig. 1. This effective energy release is 
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constrained to be smaller than 10 4 |[L3|| . Similar to Eq. ( |23| ) this leads to the lower limits 

1.2 x 10 5 for monopole annihilation 
> I 1.5 x 10 5 for ordinary strings \ x f b ( — — ) (25) 



"diV 7 


( 2-Eobs \ 


dx 


V mi / 



1.1 x 10 10 for the OTW scenario 



0.75 



m x 



A 7 (£ obs )\ 1 / jHECR(-Eobs) • GeVcm 2 secsr N 



V10 16 GeV/ V lOMpc / \ 4 x lO" 31 / 

These constraints are less stringent than the constraints Eq. (^3|) from cascade nucleosyn- 
thesis. For the OTW scenario effective power law fragmentation functions are inconsistent 
with CMBR distortions for f b > 0.11. 

It should be noted that in the superconducting string scenario there is an additional con- 
tribution to the CMBR distortions even if HECR are not produced at all. This contribution 
comes from the Sunyaev-Zeldovich effect caused by the hot gas produced around the string 
by emission of electromagnetic radiation before it reaches saturation length and potentially 
starts to emit HECR. This was discussed in Ref. [0. Our restriction to distortions caused 
by HECR alone therefore renders our constraints conservative. 



3.4 Limits from the 7-ray Background 

Electromagnetic cascades which are started at relatively low redshifts z produce an isotropic 
7-radiation in the observable energy range. The flux in this radiation puts an upper limit 
on the possible flux of ultrahigh energy particles. The most stringent constraint comes from 
the upper limit to the observed isotropic flux at E 1 ~ 200 MeV, which was reported to be 
7x lO^MeVcn^secsr)- 1 [gl] . 

The limits derived below crucially depend on the assumptions about fragmentation of 
X-particles into the usual particles like protons, pions, photons, electrons etc., and on the 
assumption about cosmological evolution of X-particle production [see Eq. (TO)]. 

We shall assume that the fragmentation function for the decay of X-particles with mass 
mx into particles i {i =p,7,e) has the form 

2 dN dN 1 ft)R . 

V x ) = -JlT\ E h ™>x) = Al -=r, (26) 



mx dx dEi \mx J E, 

where Ei is the energy of particle i and Ai is a normalization constant. For q we shall focus 
on the values between q = 1 inspired by scaling distribution in inelastic pp-scattering and 
q = 1.32 according to QCD calculations p. 

As far as evolution is concerned we shall consider two cases: (i) absence of evolution and 
(ii) the "weak" evolution, as given by Eq. (p!9|)and inspired by the development of a network 
of cosmic string loops p2|. The strong evolution with p < 1 [see Eq. (IT])] results in more 



stringent limits and we shall skip it in this paper. 

Let us first turn to the non-evolution case (i). Let the HECR flux observed at _E b s = 
2 x 10 20 eV, jHECR(-Eobs) — 4 x 10~ 31 ( GeV cm 2 sec sr) _1 , be caused by protons or 7-rays. The 
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generation function for these particles in GeV 1 cm 3 sec 1 can then be found as 

47T 

HKbs) = , fJP v jHECR^obs) , (27) 

which also leads to Eq. (|20"D . This can be extrapolated to other energies by using the fragmen- 
tation function Eq. fl26|) . In Eq. fl27j) i =p or 7, and Xi(E ohs ) is again the attenuation length 
for these particles in the CMBR field. From Eqs. ( p6|) and (|2^) one can then find the total 
energy production in form of protons, 7-rays and electrons (i =p,7,e) in GeV cm -3 sec -1 . 

The energy released in electrons and 7-rays (produced directly or through the decay of 
other particles) goes into electromagnetic cascades (the cascade energy production due to 
protons is considerably less). Using the usual quark counting one can estimate that about 
10% of the total energy release goes into electrons and thus into the cascades. The flux of 
the cascade photons can then be found as 



•cas/p \ _ c (2/3)iJ 1 q cas rp-3/2 / 9S \ 

J 7 l-frobs) - — : ; ; — 572 E , (28) 

7 4tt [ 2 + \n(E a /E x )]El' 2 7 

where E a and E x are characteristic cascade energies which for z = are given by E a ~ 
8 x 10 4 GeV and E x ~ 5.1 x 10 3 GeV, and g cas is equal to the energy release in the form of 
electrons and 7-rays. 

From Eqs. fl27|) and (§§ we find the cascade flux at £ 7 ~ 200 MeV to be 8 x 10~ 8 , 
3 x 10 -8 , and 9 x 10~ 9 ( MeV cm 2 sec sr) _1 for q = 1.1, 1.2 and 1.32, respectively, assuming 
mx = 10 16 GeV. These numbers should be compared with the observational upper limit 
7 x 10 _8 ( MeV cm 2 sec sr) -1 [61]. For q = 1.32 the predicted flux is one order of magnitude 



less. 

Let us now go over to the case of evolution (ii). The cascade limit becomes more stringent 
in this case because the cosmological epochs with large z give no contribution to the presently 
observed HECR flux at E ~ 10 20 eV, while they contribute strongly to the cascade energy 
density due to the enhanced energy release at earlier times. We shall restrict ourselves to 
the case of weak evolution here where integration over redshifts results only in a logarithmic 
factor. 

It is easy to understand the existence of a "critical" epoch (with redshift z c ) in our 
problem. It is defined as x (1 + z c ) = E x (z c ), where E 7 is a photon energy at z = and 
E x (z c ) is the turn-over energy of the cascade spectrum at redshift z c . For E^ ~ 200 MeV 
one finds z c ~ 100. If we integrate the evolution function Eq. (^) over the redshift interval 
between z = and z = z c we obtain 



h 1 ° bsJ 4tt 2 + HE a (z c )/E x (z c )) [EMY /2 " ' [ ) 

where g cas is found with the help of Eq. ( P?D and the fragmentation function Eq. ( PB|) using 
the energy transfer into p/y and e at large redshifts. 



For q = 1.1, 1.2, and 1.32, the flux Eq. (||) at £ 7 ~ 200 MeV is numerically 8 x 10 



-6 



3 x 10 6 , and 9 x 10 7 ( MeV cm 2 sec sr) 1 , respectively. For X-particle masses different from 
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mx = 10 16 GeV these fluxes have to be multiplied by (mx/10 16 GeV) 2-9 . These numbers 
are considerably higher than the upper limit 7 x 10~ 8 ( MeV cm 2 sec sr) _1 as long as mx is 
not much smaller than 10 16 GeV. These considerations can be translated into the lower limit 
q > 1.6 for the index of an assumed power law injection. 

Note that Chi et al. [63] derived similar limits by considering cascade development in the 
CMBR and in the infrared and starlight fields. These limits depend to some extent on the 
history and intensity of these less well known backgrounds. However, in the case of "weak 
evolution" of TDs considered here the comparatively strong injection at high redshifts leads 
to cascading probably mostly in the CMBR, whereas the authors of Ref. |64| were more 



concerned with low redshift injection where these other backgrounds are more important. 

As a conclusion we claim that for a fragmentation function of the form of Eq. (|26"D with 
reasonable values for q, 1 < q < 1.32, the explanation of observed HECR at E > 10 20 eV as 
protons or 7-rays from the decay of GUT scale X-particles with mx — ^gut — 10 16 GeV 
is incompatible even with the "weak" cosmological evolution of their production. The non- 
evolution case is not severely constrained by these arguments. 



4 Conclusions 

We have discussed limits on cosmic high energy particle injection derived from 4 He photo- 
disintegration, CMBR distortions and the diffuse 7-ray background. We have found that the 
nucleosynthesis limits give the most stringent constraints for epochs with redshift z > 5 x 10 3 
whereas at lower redshifts particle injection is predominantly limited by its contribution to 
the diffuse 7-ray background (see Fig. 1). These constraints were applied to topological 
defects potentially radiating supermassive GUT scale ("X") particles which subsequently 
decay into high energy leptons and hadrons. The history of high energy particle injection 
is more or less determined within these defect models. The model dependent parameters to 
be fixed are the number density of X-particles radiated within unit time and the effective 
fragmentation function for the decay products of these X-particles. We have assumed that 
the flux of these decay products contributes significantly to the present day observed HECR 
flux. This allowed us to formulate our constraints as lower limits on the fractional energy re- 
lease at HECR energies (~ 10 20 eV) which is mainly determined by the 7-ray fragmentation 
function. We have found that for reasonable 7-ray fragmentation functions superconducting 
strings can not explain the HECR flux without violating at least the bound coming from 
4 He-photodisintegration. In contrast, magnetic monopole and ordinary cosmic string models 
producing observable HECR fluxes are most severely constrained, but not yet ruled out, by 
their contribution to the diffuse 7-ray background. 

In the second part of the paper we have studied the possibility that the presently ob- 
served deuterium has been produced by an epoch of 4 He-photodisintegration subsequent to 
a standard nucleosynthesis scenario. Such an epoch may have been initiated by the decay of 
particles, the annihilation of topological defects, or, in general, the production of energetic 
7-rays by any source. We have found that only a small fraction (< 10%) of the observed 
deuterium may have its origin in the process of 4 He-photodisintegration since, otherwise, 
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anomalously large primordial ( 3 He/ 2 H)-ratios would result. A larger fraction of the pri- 
mordial deuterium contributed by this process would require either standard assumptions of 
chemical evolution to break down or the existence of 7-ray sources in the early universe which 
radiate with extremely "soft" 7-ray energy spectra. We have shown that a scenario which 
employs massive black holes to reprocess the light element abundances from a standard big 
bang nucleosynthesis process [[Uj is in conflict with 2 H and 3 He observations. We have also 
used the anomaly in the ( 3 He/ 2 H)-ratios produced during 4 He-photodisintegration to slightly 
tighten constraints on the abundances and parameters of decaying particles and topological 
defects. 
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Figure Captions 

Figure 1: Maximal energy release in units of the CMBR energy density allowed by the 
constraints from the observed 7-ray background at 200MeV (dotted curve), CMBR distor- 
tions (dashed curve, from Ref. |13[), and 4 He-photodisintegration as a function of redshift z. 
These bounds apply for instantaneous energy release at the specified redshift epoch. The 
logarithm is to the basis 10. 
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